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Dark Matter annihilation or decay can affect the anisotropy of the cosmic microwave background 
(CMB). Therefore, the CMB data can be used to constrain the properties of a dark matter particle. 

In this work, we use the new CMB data obtained by the Planck satellite to investigate the limits 
on the basic parameters of a dark matter particle. The parameters are the dark matter mass (m^) 
and the thermally averaged cross section ((ctu)) for dark matter annihilation and the decay rate 
(T) (or lifetime r = l/T) for dark matter decay. For dark matter annihilation we also consider the 
impact of the structure formation process which is neglected by the recent work. We find that for DM 
annihilation, the constraints on the parameters are fann = {cm) < 0.16x 10“^®cm®s“^GeV“^(or 
fann < 0.89 X 10~®m®s“^kg“^, 95% C.L.). For DM decay, the constraints on the decay rate are 
r < 0.28 X 10-2®s-1(95% C.L.). 


I. INTRODUCTION 


Dark matter (DM), the main components of the Uni¬ 
verse, confirmed by many observations still keeps its 
nature mysterious M. Its most widely accepted 
model, the weakly interacting massive particles (WIMPs) 
model, claims DM particles can self-annihilate into stan¬ 
dard model particles, such as photons, electrons, and 
positrons [I| , which might be observed by various exper¬ 
iments, such as PAMALA and AMS-2 Q. During the 
evolution of the Universe, these particles produced by 
DM annihilation interact with the medium of the Uni¬ 
verse For example, the photons produced through 

DM annihilation can ionize the hydrogens formed in the 
epoch of the recombination before appearance of the first 
stars. The changes of the ionization will be reflected 
in the anisotropy of the cosmic microwave background 
(CMB) (7l-[Tlj|. Therefore, the observation data of CMB 
can be used to investigate the nature of DM particles. 
Recently, the authors of [ll| used the Planck data to in¬ 
vestigate the limits on the DM parameters, the mass m^, 
and the thermally averaged cross section {av). In that 
paper, they considered the uniform distribution of DM 
while neglecting the structure formation process of the 
Universe which claimed that the DM halos were formed 
in the redshift z ^ 100. Moreover, the subhalos or sub¬ 
subhalos are also formed in DM halos [l2|. The DM an¬ 
nihilation rate can be enhanced in all these DM halos. In 
this work, we consider these effects. In addition to anni¬ 
hilation, DM particles are not stable in some models, and 
they can decay into the standard model particles [l^ [l^ . 
In this paper, we also use the Planck data to get the con¬ 
straints on the decay rate (or lifetime) of DM. Not setting 
any specific DM model, the results of this work can be 
applied widely to those DM particles can annihilate or 
decay. 

This paper is organized as follows: In Sec. II, we briefly 
show how the DM particles’ annihilation and decay affect 
the evolution of the Universe and we give the constraints 
on the basic parameters of DM by use of the Planck data. 
We conclude in Sec. III. 


II. THE IMPACT OF DM ON THE EVOLUTION 
OF THE UNIVERSE AND CONSTRAINTS ON 
THE BASIC PARAMETERS 

DM particles can affect the Universe through the inter¬ 
action between the medium of the Universe and the pro¬ 
ductions of DM annihilation or decay. The two main ef¬ 
fects are to ionize the hydrogens and heat the medium . 
The changes of ionization with the redshift are governed 
by the equation @,0 


{lPz)-^ = [Rs{z)-I,{z)-I^{z)\ (1) 

where is the ionization rate due to the DM, and 

Rs{z) and Is{z) are the standard recombination rate and 
ionization rate by the standard sources, respectively. The 
ionization rate by the DM can be be written as 




Xif {z)^ann 


UbEb 


for the DM annihilation and 


( 2 ) 




Xif {^)^dec 


rub 

Ulb ribEb 


( 3 ) 


for the DJVI decay. Pann and Pdec (or Puec — 'i~dec 

is the lifetime ) are the DM annihilation and decay rate, 
Ub is the number density of baryon, and Eb = 13.6eV is 
the ionization energy. 

If the structure formation effect is included, the DM 
annihilation rate can be written as 


=fannPcri^ci^ + z)^B{z) (4) 

where B{z) is the ’boost factor’ (or ’clumping factor’) 
due to the structure formation effect; for more detailed 
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discussions, one can see Ref. [l^. pcri is the criti¬ 
cal density of the Universe, farm is a new parameter, 
which is the combination of the basic parameters of DM, 

farm — 

For the boost factor B{z), we followed the Ref.(l5l| and 
used the form as 

where dn/dM is the mass function of DM halos and 
we use the Press-Schechter formalism for our calcula¬ 
tions [l^. For the DM halos, we use the NFW (Navarro- 
Frenk-White) density profile. It is also found that there 
are many sub- and sub-substructures in DM halos [l^ . 
The DM annihilation rate can be enhanced in these sub¬ 
halos. In this paper, we include these subhalos while 
neglecting the contributions from the sub-subhalos. We 
set the smallest mass of DM halo as ^ IO^^Mq [l^. We 
adopt that 10% mass of DM halos is in the form of sub¬ 
halos. We use the power law form of the mass function 
~ for the subhalos [l^. The total boost factor of 

DM halos including the subhalos is [Tsi l 


structure formation. In this work, we include this effect, 
and we also get the constraints on the decay rate (or 
lifetime) for the DM decay. Because the limits are inde¬ 
pendent on any specific DM models, the results can be 
used widely for many DM annihilation or decay models. 

We modified the public code RECFAST 0to included 
the DM effect, and we used the public code CosmoMC0 
to get the constraints on the parameters. We consider 
six cosmological parameters and a new parameter 

,Dc/i ,d,T,TZs, Ag,(8) 

for the DM annihilation and 

{nbh'^,nch'^,9,T,ns,As,Fdec}, (9) 

for the DM decay. Here are the new parameters: Fann = 
f{z)fann = f{z){av)/my^ and Fdec = Tdecfiz). 

The final constraints on the parameters are given in 
Table [T] for the DM annihilation and Table HIl for the DM 
decay. For making comparisons, the results for the case 
of smooth DM distribution are given too. In this case, 
the ’boost factor’ is B{z) = 1. From these results, it can 
be seen that for f{z) = I, the 95% upper limits are 


Btotal — 1 T {Bhalos I) T {Bgabhalos 1) (6) 

In equations ([2]) and Xi stands for the fraction of 
the energy which contributes to the ionization and is dis¬ 
cussed first in Ref. [l^. Here we adopt the form given 
by Ref. Xi = (1 ~ Xe)/i- It should be noticed that 
although this form has been used frequently in pre vious 
works Hi, 13 , it is not accurate. In the Ref. |l9l |. the 
difference of the constraints on DM parameters occured 
both in the approximate and the accurate have been in¬ 
vestigated and only slightly differences of the upper limits 
for the present null detection of DM found. The similar 
discussions are also present in the Ref. [l^. f{z) is the 
fraction of the energy which deposits to the medium of 
the Universe. It is different for different annihilation or 
decay channels and is a function of the redshift. Here we 
treat it as a free parameter. For more detailed discus¬ 
sions one can see Ref. [TTj | . In addition to the ionization, 
another important effect of DM to the evolution of the 
Universe is to heat the medium, e.g, the baryonic gas. 
The evolution of the gas temperature can be written as 




SaraRT^rnb 

'imeCH (z) I -I- /i/e + Xe 


{Tb — Tcmb) 


2 _ 

3kBH{z) 1 -I- fne + Xe 


+ 2Tb 


(7) 


where is similar to and stands for the heating rate 
to the medium by the DM. Following Ref. Q, we use the 
form as = (1 + 2xe)/3. 

Recently, the authors of m used the Planck data to 
get the constraints on the DM parameters for the anni¬ 
hilation case, but they did not consider the effect of the 


< 0.16(0.24) X lO-^^cm^s-^GeV^ (10) 

for the DM annihilation case (0.24 is the case of smooth 
distribution) and 


F < 0.28 X 10“®^s“^ (11) 

for the DM decay case. 

For the DM annihilation case, our results are consis¬ 
tent with Ref. [ni(e.g., Table H). For this point, we can 
convert our results as 


— < 0.16(0.24) X lO-^Scm^s-^GeV-i 

= 0.89(1.34) X lO’V^s-^kg-i (12) 

We also plot the ID and 2D probability distributions 
in Figs. [T] and [2] for the DM annihilation (including the 
smooth distribution case, dotted lines) and Figs. |3] and 
m for the DM decay, respectively. From these plots, one 
can find that the correlation between the Farm parameter 
and the cosmological parameters is stronger for the DM 
annihilation than that of F^ec for the DM decay. The 
main reason is that the DM annihilation rate is propor¬ 
tional to the number density square, so the effects due to 


^ http: //camb.info 
^ http://cosmologist.info/cosmomc/ 
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the annihilation are very strong during the recombina¬ 
tion. For the DM annihilation, the differences between 
the clumpy and smooth DM distribution cases are not 
so huge. These results indicate once again that the lim¬ 
its on the DM parameters are mainly from the epoch of 
recombination. 


III. SUMMARY AND DISCUSSION 

In this work, we used the new data from the Planck 
satellite to investigate the limits on the DM basic pa¬ 
rameters for the annihilation and decay. By consider¬ 
ing the structure formation effect for the DM annihi¬ 
lation case, we found that the constraints on the fann 
parameter are fann < 0.16(0.24) x 10“^®cm^s“^GeV“^ 
or fann < 0.89(1.34) X 10-6m3s-ikg-i(95% C.L.). For 
the DM decay, the constraints on the decay rate are 
F < 0.28 X 10-25s-1(95% C.L.). 

As mentioned in the Sec.II, for the clumpy DM distri¬ 
bution, the smallest mass of DM halo is set as ^ 10~^Mq, 
which is different for different DM models. In theory for 
WIMPs DM, this value ranges from to 10~^Mq 

for typical kinetic decoupling temperatures. From the re¬ 
sults of current numerical simulations, the typical small¬ 
est mass of DM halos is ^ 10 ®Mq. In Ref. [T^, the 
authors discussed the effects on the ’boost factor’ for the 
different values of the smallest DM halos. They found 
that there are differences of ^ 2 orders of magnitude 
of the ’boost factor’ for DM halo mass and 

10“"^Mq at z ~ 50 (upper panel of Fig. 1 in Ref. [l5l|b 
For DM halo mass 10 ®Mq, the differences of ~ 5 orders 
of magnitude differences are present at z ~ 20 compared 
with the DM smooth distribution. Therefore, the largest 
differences usually appear in nearby universe, and it is 


believed that the changes of limits on the DM param¬ 
eters are slight if one change the values of the smallest 
mass of DM halos. 

Another factor which can affect the limits on the DM 
is the density profile of DM halos. In this work, we have 
used the NFW profile, which is well in fitting many ob¬ 
servations data. In addition there are still many other 
observations or N-body simulations which are favored by 
the other profiles, such as Einasto profile [2ll - l^ . which 
are slightly different from that of NFW profile for the 
final constrains. 

One point in this work that should be noticed is that 
we have set f{z) as a free parameter, and for the final 
constraints we have set /(z) = 1, which means that all 
the energy produced by the DM annihilation or decay has 
deposited into the medium of the Universe. In Ref. 0, 
the dependence of /(z) on the redshift and different an¬ 
nihilation channels were discussed by the authors. It can 
be seen that the final constraints are slightly different 
(Table II of Ref. 0. 

NOTE: the latest eonstraints on the DM parameters 
for annihilation can be found in the paper of Planck Col¬ 
laboration (arXiv:1502.01589). 
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FIG. 1. The marginalized probability distribution function of parameters for the DM annihilation case. The solid and dotted 
lines correspond to the clumpy and smooth DM distribution respectively. 
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TABLE 1. Posterior constraints on the parameters for the DM annihilation. The first and second lines of every item correspond 
to the clumpy and smooth DM distribution respectively. 


Parameter 



1000 

T 

Us 

ln{lG^°As) 

^^^^(lO-^^m^s-iGeV-i) 

Mean 

0.02207 

0.1196 

1.0412 

0.088 

0.962 

3.09 

0.053 


0.02209 

0.1195 

1.0413 

0.089 

0.963 

3.10 

0.070 

2a low 

0.02155 

0.1147 

1.0401 

0.061 

0.948 

3.05 

0 


0.02152 

0.1140 

1.0401 

0.062 

0.949 

3.05 

0 

2a up 

0.02260 

0.1251 

1.0425 

0.1149 

0.977 

3.15 

0.16 


0.02265 

0.1244 

1.0425 

0.1154 

0.978 

3.15 

0.24 



TABLE 11. Posterior constraints 

on the parameters for the DM decay. 


Parameter 


Q.ch?' 

1000 

r 

Us 

ln{lQ^°As) 

Edec(10- = ®S-1) 

Mean 

0.02205 

0.1199 

1.0413 

0.088 

0.960 

3.09 

0.077 

2a low 

0.02152 

0.1147 

1.0401 

0.064 

0.94 

3.04 

0 

2a up 

0.02262 

0.1251 

1.0425 

0.115 

0.97 

3.14 

0.28 








FIG. 2. The 2D contours distribution function of parameters for the DM annihilation case(68% and 95% confidence level). The 
solid (black) and dotted (red) lines correspond to the clumpy and smooth DM distribution respectively. 
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FIG. 3. The marginalized probability distribution function of parameters for the DM decay case. 
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FIG. 4. The 2D contours distribution function of parameters for the DM annihilation case(68%, 95%, and 99% confidence 
level). 















